Short to long-range charge-transfer excitations in the 
zincbacteriochlorin-bacteriochlorin complex: a Bethe-Salpeter study 
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We study using the Bethe-Salpeter formalism the excitation energies of the zincbacteriochlorin- 
bacteriochlorin dyad, a paradigmatic photosynthetic complex. In great contrast with standard time- 
dependent density functional theory calculations with (semi)local kernels, charge transfer excitations 
are correctly located above the intramolecular Q-bands transitions found to be in excellent agreement 
with experiment. Further, the asymptotic Coulomb behavior towards the true quasiparticle gap for 
charge transfer excitations at long distance is correctly reproduced, showing that the present scheme 
allows to study with the same accuracy intramolecular and charge transfer excitations at various 
spatial range and screening environment without any adjustable parameter. 

PACS numbers: 71.15.Qe,78.67.-n,78.40.Mc 
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Photoinduccd charge transfer excitations, namely the 
jump upon photon absorption of an electron from a donor 
to an acceptor site, is a fundamental process that gov- 
erns photosynthetic processes in plants and bacteriasji 
or the quantum efficiency in organic or hybrid photo- 
voltaic cells.— Such non-local excitations are also an im- 
portant current theoretical issue since it is now well rec- 
ognized that the time-dependent density functional the- 
ory (TDDFT)^ encounters severe problems to describe 
such excitations when standard (semi)local kernels, or 
even hybrid kernels mixing some amount of exact ex- 
change, are being used^ Besides organic systems, similar 
problems have been identified in the case of extended 
Wannicr excitons in semiconductors where the large ef- 
fective excitonic radius leads to a weak average overlap 
between the hole and the electron^ 

The bactcriochlorin molecule is closely related to the 
magnesium-containing bacteriochlorophyll system. Due 
to its importance as a paradigmatic photosynthetic com- 
plex, and as one of the earliest charge-transfer system 
for which the TDDFT difficulties have been unraveled 
and discussed,— the zinc-bacteriochlorin/bacteriochlorin 
(ZnBC-BC) complex (see Fig. [T]) has been studied by a 
variety of approaches, including TDDFT with local^^ 
hybrid^ Coulomb attenuated hybrid 7 - functionals, con- 
strained ASCF DFT calculations^ and quantum chem- 
istry many-body wavefunctions techniques such as a com- 
bination of ASCF DFT and single excitation configura- 
tion interaction (CIS) technique^ or a more elaborate 
CIS(D) approach including various scaled perturbativc 
double-excitation correlation corrections^ 

In two recent studies ; 10 ' 11 charge transfer (CT) ex- 
citations in small donor/acceptor complexes, combining 
tctracyanocthylcnc with acene derivatives, were studied 
with the GW approximation and Bethe-Salpeter (BSE) 
equation within many-body perturbation theory.— Ex- 
cellent agreement with gas phase experiments^ was ob- 
tained for the lowest CT excitation energy with a mean 
absolute error of about 0.1 eV.— Such an accuracy com- 
pared well with recent TDDFT calculations with opti- 



mized range-separated functionals^ while TDDFT cal- 
culations with standard PBE or even non-local B3LYP 
kernel a 15 ' 16 were shown to lead to discrepancies of several 
eV with CT states located at much too low energy*^ 

In this work, we study the optical absorption spectrum 
of the ZnBC-BC complex using the GW-BSE many- 
body perturbation theories. We show that intramolec- 
ular excitations arc in excellent agreement with experi- 
ment, and that the charge-transfer excitations are cor- 
rectly located above the intramolecular Q-bands transi- 
tions, in great contrast with TDLDA calculations which 
locate the charge-transfer excitations about 1.7 eV below 
the GW-BSE value. It is further shown that the GW- 
BSE framework correctly reproduces the long-range en- 
ergy behavior of the charge transfer excitations, a feature 
shared only with the scaled CIS(D) correlated quantum 
chemistry approach. The importance of going beyond 
the Tamm-Dancov approximation is further discussed. 

Our calculations arc based on a recently developed 
auxiliary even-tempered Gaussian-basis^ implementa- 
tion of the GW^r^- and BSE22r— formalisms, the Fi- 
esta package, implementing contour deformation tech- 
niques for the correlation contribution to the self- 
energy. 11 ' 27 The needed two-body operators such as 
the susceptibilities, bare/screened Coulomb potentials 
and self-energies, are expressed on an atom-centered 
Gaussian-basis containing six exp(—ar 2 ) gaussian for the 
radial part of each (s, p, d)-channcls, with an even tem- 
pered distribution of localization coefficient a ranging 
from 0.1 to 3.2 a.u. Such a basis derives from previous 
studies^ 7 - - — but with additional diffuse orbitals. 

The needed starting single-particle states are taken to 
be the Kohn-Sham DFT/LDA cigenstates generated by 
the Siesta package^ , using a large triple-zcta plus double 
polarization basis (TZDP)^i It was shown recently that 
the combination of a TZDP basis for describing Kohn- 
Sham occupied and virtual states, with the above de- 
scribed auxiliary bases, was leading to charge-transfer ex- 
citation energies in excellent agreement with planewave 
GW-BSE calculations^ For the isolated ZnBC and BC 
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FIG. 1: (Color online) Symbolic representation of 

the ( 1 ,4)-phenylene-linked zincbacteriochlorin-bacteriochlorin 
complex. Isocontour representation of the (grey) hole- 
averaged electron and (blue) electron-average hole distribu- 
tions (see text) for the lowest (a) ZnBC Q x excitation and 
(b) ZnBC — > BC charge-transfer excitation. 



molecules, for which experimental data are available ; 32 ' 33 
we relax the molecules at the all-electron B3LYP 6- 
311G(d) levels However, for the ZnBC-BC bridged com- 
plex, we adopt the geometry provided in Ref. |H, and used 
in subsequent studies* 7 - - — as a mean to properly compare 
to other computational approaches. 

In the present many-body framework, one first ob- 
tains accurate quasiparticle energies (occupied and vir- 
tual electronic energy levels) by using the GW formal- 
ism. In the case of gas phase organic molecules such 
as fullerenes, pentacene, PTCDA or porphyrins, it was 
shown in several recent studie d 27 ' 35- — that the stan- 
dard single-shot Go Wo calculations based on input Kohn- 
Sham eigenstates with standard local functionals tends to 
yield too small HOMO-LUMO gaps, even though already 
much better than the starting Kohn-Sham values. Better 
results could be consistently obtained for the quasiparti- 
cle spectrum of small molecules either by starting from 
Hartree-Foc k 27 ' 38 " — , or hybrid functionals 41 - eigenstates, 
or by performing a full self-consistency on the quasipar- 
ticle eigenstates , 39 ' 40 or a limited self-consistency on the 
eigenvalues only ; 11 ' 27 ' 29 ' 37 ' 42 allowing to remove the de- 
pendency of the final quasiparticle energies on the start- 
ing point eigenstates. This is such an approach that 
we adopt in the present study, showing here below as 
in Ref. [ll| that excellent agreement with experiment can 
be obtained for the excitation energies. 

In an optical absorption spectrum, the quasiparticle 
gap is reduced by the electron-hole interaction. In the 
present MBPT formalism, this excitonic interaction is 
accounted for by the BSE equations. Namely, the neu- 
tral excitation energies can be obtained as the eigen- 



values of the BSE H e h Hamiltonian^ 2 - — expressed in 
the 4>i(r)4>j(r') product basis of the unoccupied 4>f and 
occupied 4>j single-particle states. Of importance for 
the upcoming discussion, we write the so-called direct 
term involving the (statically) screened Coulomb poten- 
tial W(r, r'), noticing that hole and electron states are 
not taken at the same space position: 

=~j drdr>f(r)^(r')W(r,r')^(r)^(r'), 

As compiled in Table Q] we first verify that our GW-BSE 
low-lying Q x and Q y transition^ 3 - energies for the isolated 
monomers fall within 0.1 eV of the available experimen- 
tal values. Consistently with previous observations ; 28 ' 44 
we find that the diagonalization of the full BSE Hamil- 
tonian leads to a red shift which can be as large as 0.35 
cV as compared to the GW-BSE results in the Tamm- 
Damcov (TDA) approximation, bringing the calculated 
transitions in excellent agreement with experiment. Such 
a good agreement with available experimental data can 
be taken as an indication of the reliability of the present 
formalism and of the specific implementation aspects. 

Concerning the TDLDA calculations, performed with 
the same Kohn-Sham states and auxiliary basis, we ob- 
serve that the agreement with experiment is also satis- 
factory, even though not as good as the BSE values. Our 
calculated TDLDA results compare very well with for- 
mer TD-BLYP calculations.— i 7 - Again, this allows to ver- 
ify that the pseudopotential approximation and the bases 
we use do not lead to significant errors as compared to all- 
electron calculations with standard quantum chemistry 
basis. For intramolecular transitions with large overlap 
between hole and electron states, the TDDFT approach 
with standard kernels is a reliable framework. 

We now come to the central point of this study, namely 
the ZnBC-BC CT excitation energies. The charge- 
transfer nature of a given transition can be easily identi- 
fied by analyzing the weight of the two-body BSE eigen- 
states on the 4>fumo an d ^homo one-body eigenstates lo- 
calized on either molecule. The expectation value of the 
electron/hole density operator <5(r — r e /r^,) on a given 
two-body BSE eigenstate allows further to build the cor- 
responding hole/electron-averaged electron/hole distri- 
bution. Such a representation is provided in Fig. [ljb) 
for the ZnBC— >BC charge-transfer state obtained at the 
GW-BSE (full) level. The charge-transfer nature of such 
an excited state is apparent. 

As shown in Table HI TD-LDA severely underestimates 
the CT excitation energies which are found to lie signifi- 
cantly below the intramolecular Q-bands. Identical con- 
clusions were found with the BLYP functional.— i 7 - Such 
a behavior can be understood from the analysis of the 
TDDFT electron-hole coupling terms which vanish for 
transitions between non-overlapping electron and hole 
states in the case of (semi)local exchange-correlation ker- 
nel. One is then left with the non-interacting diagonal 
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TD-DFT (TDA) 


GW-BSE 


Exp. c ' d 


Transitions 


BLYP a/i> 


LDA B3LYP b 


CAM-B3LYP 6 


TDA full 










Isolated BC and ZnBC monomers 




7T-7T* Q 


> x ZnBC 


2.07/2.04 


2.09 


2.05 


1.87 


1.94 1.59 


1.65 


7T-7T* Q 


>« BC 


2.10/2.10 


2.10 


2.12 


1.92 


1.90 1.63 


1.6 


7T-7T* Q 


)y BC 


2.39/2.40 


2.45 


2.54 


2.53 


2.38 2.24 


2.3 


7T-7T* Q 


> H ZnBC 


2.44/2.43 


2.47 


2.60 


2.59 


2.29 2.27 


2.2 








ZnBC-BC phenylene-bridged complex 




7T-7T* Q 


ZnBC 


2.05/2.05 


2.08 


2.41 


1.89 


1.99 1.69 




7T-7T* (J. 


>« BC 


2.09/2.09 


2.10 


2.43 


1.93 


2.04 1.73 




7T-7T* Q 


»» BC 


2.38/2.37 


2.46 


2.79 


2.49 


2.47 2.35 




7T-7T* (J. 


>„ ZnBC 


2.42/2.42 


2.44 


2.83 


2.54 


2.49 2.34 




CT ZnBC -> BC 


1.33/1.33 


1.25 


1.96 


2.87 


2.95 2.95 




CT BC ->■ ZnBC 


1.46/1.46 


1.39 


2.12 


3.04 


3.13 3.13 





TABLE I: Calculated singlet transition energies (in eV) for the isolated monomers and the (l,4)-phenylene-linked ZnBC-BC 
complex. The present TD-LDA (TDA), GW-BSE (TDA) and GW-BSE (full) results are compared to previous TD-DFT 
calculations with the BLYP and B3LYP functionals. (TDA) means Tamm-Dancov approximation, while (full) means full 
diagonalization mixing resonant and antiresonant transitions. a Ref. |. 6 Ref.Q. c ' d Refs. \MM- 



part, namely the too small Kohn-Sham HOMO-LUMO 
gap. Such a situation is slightly improved with the 
B3LYP functional due to its 20% of exact exchange, but 
with a residual error as large as 1 cV. This cancellation 
of the interacting term for non-local charge-transfer exci- 
tations does not occur in the BSE formalism since in the 
HfYki* interaction term, hole and electron states are not 
taken at the same position in space and are connected by 
the non-local screened Coulomb potential W(r, r'). 

The GW-BSE CT excitation energies lie well above the 
intramolecular Q-bands, consistently with TDDFT cal- 
culations with optimized (parametrized) CAM-B3LYP 
functional^ and CIS(D) calculations^ (see analysis be- 
low). Concerning the CAM-B3LYP results, with CT 
states located less than 0.1 eV below the GW-BSE 
one, to be compared to the ~ 1.7 eV discrepancy with 
TDLDA, the agreement is very good considering that 
the needed (a, f3) scaling parameters 4 ^ have been trained 
on a very different set of molecules, the so-called G2 
set, and on very different properties, namely ionization 
energies. These results confirm the agreement already 
found between GW-BSE and a TDDFT calculation with 
another optimized range-separated functional (the BNL 
functional 4 ^) in the study of small TCNE-accncs com- 
plexes, both calculations coming in close agreement with 
experiment for CT states . 11 ' 14 

The ability of the BSE approach to describe charge- 
transfer excitations can be further illustrated by study- 
ing the long range limit where the BC and ZnBC 
units are separated by removing the phenylene bridge. 
The distance R between the two monomers in the un- 
bridged model dimcr is defined in Fig. 2a (Inset). R <~ 
5.84 A is the corresponding distance between the two 
monomers in the phenylene-bridged complex. Our re- 
sults are represented in Fig. O As expected, intramolec- 
ular charge transfer excitations are independent of the 
distance between the molecules. On the contrary, in 
the large distance limit, the CT exciton binding ener- 
gies are found to scale like 1/D, where D is the dis- 



tance between the R\ and i?2 molecule centres, as in- 
dicated in Fig. [2b (thick grey line). In the large D limit, 
the direct interaction terms converge to: Hf? ^ ct ~ 

-W (Ri, TL 2 )Ql k Qjii where W(Ri,R 2 ) reduces to the 
bare 1 / 1 Ri — R2 1 Coulomb potential for the two molecules 
in the vaccuumil and Qj t =< <f>\\4>\. >■ With Q e ik = Sik 
and Qjj = Sji, there is no mixing with higher subbands, 
so that the lowest CT state is primarily composed of the 
HOMO and LUMO eigenstates in the long distance limit, 
as verified by analyzing the BSE eigenstates. 

Such (1/D) behavior of the CT excitonic binding en- 
ergy is the correct asymptotic limit of an electrostatic 
interaction between opposite charges. In this limit, the 
CT excitation energies correctly converge towards the 
non-interacting (diagonal) part of the BSE Hamiltonian, 
namely the true quasiparticlc gap between the concerned 
hole and electron states as given by the GW calcula- 
tion. The ability of the GW-BSE approach to correctly 
describe without adjustable parameter both intramolecu- 
lar and short to long-range charge transfer excitations 
energies, both infinite (molecular) and extended (infinite 
solids) systems^ is a important feature as discussed 
here below by comparison with previous studies. 

As emphasized above, TD-DFT calculations with lo- 
cal kernels cannot reproduce the long-range electron- 
hole interaction: the CT energies remain constant and 
equal to the Kohn-Sham band gap. This is now well 
documented^&i and we clearly verify that point with our 
TDLDA data (see dotted lines below 1.5 cV in Fig. 
As a dramatic improvement over local kernels, the TD- 
CAM-B3LYP approach^ provides a much better agree- 
ment with GW-BSE calculations in the short-range limit 
of the bridged dimer (see above). However, as already 
noticed in Ref. 0, the (a + /?)/D Ion g-range behavior, 
with (a + /3)=0.65,— cannot reproduce the correct 1/D 
long-range limit. Similarly, the asymptotic limit of the 
TD-B3LYP approach would reduce to 0.2/D, thanks to 
the 20% of exact exchange. This discussion clearly un- 
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FIG. 2: (Color online) GW-BSE excitation energies as a 
function of the ZnBC to BC distance R in the absence of 
the phenylene bridge as indicated on the atomic representa- 
tion. The Inset indicates the intramolecular or charge-transfer 
nature of the represented GW-USE transitions. The lowest 
ZnBC — > BC CT excitation energy is compared to the {Eq^L ~ 
e 2 /D) "Mulliken" limit (thick grey line), with Eg™ the dimer 
GW quasiparticule band gap in the large D limit, where D 
is the distance between the two molecule centers. For sake 
of comparison, the TDLDA values for the two low lying CT 
excitons are indicated by the two non-dispersive dotted lines 
below 1.5 eV. The two vertical arrows indicate the difference 
between the TDLDA and GW-BSE values. For completness, 
the higher lying intramolecular Soret B-transitions and CT 
excitations, not detailed in the Table, are also reproduced. 



dcrlines the difficulty of working with fixed parameters 
functionals performing equally well at any range and in 
any screening environment (vaccuum, solvant, etc.)^ 

As a final comparison, we now analyse the difference 
with the hybrid single configuration interaction (CIS) ap- 
proach proposed in the seminal Ref. 0. In the TDA ap- 
proximation, valid for CT states in the long-range limit, 
CIS is equivalent to time-dependent Hartree-Fock show- 
ing the correct 1/D asymptotic limit for charge transfer- 
excitations, but with an incorrect asymptotic limit given 
by the donor-HOMO to acceptor-LUMO gap calculated 
in the absence of correlations. The neglect of correla- 



tion in TD-HF was corrected by a distance independent 
term obtained from DFT ASCF calculations at large dis- 
tance. For i?=5.85 A, the R distance in the true bridged 
configuration, this scheme was shown to lead a low-lying 
ZnBC -> BC CT excitation energies of 3.75 eV. This is 
~ 0.6 eV larger than our GW-BSE 3.16 eV value found 
for the model dimer at the same distance. 

As analyzed in Ref. |9|, the assumption of a distance 
independent correlation correction to CIS leads to some- 
how too large excitation energies. This problem was ad- 
dressed in a recent study^ where scaled (parametrized) 
perturbative CIS(D) double excitation corrections to CIS 
were introduced to account for correlations^ reducing 
the discrepancy with our BSE calculations to about 0.25 
eV.— While differences in basis size may possibly ex- 
plain part of the small residual discrepancy^ it remains 
that the comparison of quantum chemistry post-Hartree- 
Fock correlated techniques, such as second-order M0ller- 
Plesset (MP2) methods (see Note [49h . with the present 
GVF-BSE approach, is a current challenge, with MP2 
techniques facing scaling and divergency problems for 
extended systems with large polarizability— while the 
GVF-BSE approach on the contrary still needs futher vali- 
dation for finite size molecules. The present study clearly 
aims at contributing to that important goal. 

In conclusion, we have performed an ab initio many- 
body GVF-BSE perturbation theory analysis of the op- 
tical transition energies for the ZnBC-BC complex. For 
the isolated monomers, the calculated GVF-BSE Q-bands 
transition energies are found to be in remarkable agree- 
ment with experiment, in particular if one goes beyond 
the Tamm-Dancov approximation. In the case of the 
ZnBC-BC complex, the GIF-BSE calculations correctly 
locate the charge-transfer excitations above the monomer 
Q-bands transitions. Our CT excitation energies are 
found to be within less than 0.1 eV from parametrized 
TDDFT-CAM-B3LYP values for the short range charge 
transfer in the bridged configuration, but only the GW- 
BSE formalism can reproduce the correct long distance 
asymptotic limit. The possibility to study on the same 
footing, namely without any system-dependent parame- 
ter, charge tranfer excitation energies at various ranges, 
both in the vaccuum or in a screening environment (solid, 
solvant, etc.), opens the way to important developments 
in the study of charge-transfer excitations and energy 
transfer processes at stake in photosynthetic processes or 
in organic photovoltaic cells. 
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